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Foreword 


This  paper  was  prepared  for  presentation  at  the  Defense  Atomic  Sup¬ 
port  Agency  (DASA)  Strategic  Structures  Vulnerability/Hardening  Long  Range 
Planning  Meeting  held  at  the  U.  S.  Army  Engineer  Waterways  Experiment  Sta¬ 
tion  (WES ) ,  Vicksburg,  Mississippi,  14-16  January  1969* 

The  laboratory  tests  and  the  analysis  procedures  described  herein 
were  developed  in  conjunction  with  research  on  propagation  of  ground  shock 
through  earth  media  being  conducted  by  personnel  of  the  Soils  Division, 

WES,  for  DASA. 

This  report  was  prepared  and  presented  by  Mr.  J.  G.  Jackson,  Jr., 
Chief,  Impulse  Loads  Section,  Soil  Dynamics  Branch,  Soils  Division,  WES. 

Mr.  R.  W.  Cunny  was  Chief  of  the  Soil  Dynamics  Branch  and  Mr.  A.  A.  Maxwell 
was  Acting  Chief  of  the  Soils  Division  during  the  preparation  and  publica¬ 
tion  of  this  report.  Director  of  the  WES  was  COL  Levi  A.  Brown,  CE; 
Technical  Directors  were  Messrs.  J.  B.  Tiffany  and  F.  R.  Brown. 
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Conversion  Factors,  British  to  Metric  Units  of  Measurement 


British  units  of  measurement  used  in  this  report  can  be  converted  to  metric 


units  as  follows : 

Multiply 

By 

To  Obtain 

feet 

0.3048 

meters 

pounds  per  square  inch 

0.070307 

kilograms  per  square  centimeter 

pounds  per  cubic  foot 

It. 0185 

kilograms  per  cubic  meter 

Summary 

Computer  codes  used  in  the  solution  of  free-field  ground  shock  prob¬ 
lems  are  based  on  mathematically  defined  constitutive  models.  Quantitative 
input  for  these  constitutive  models  is  primarily  based  on  laboratory  test 
data,  but  extensive  analysis  and  numerous  assumptions  are  required  to  con¬ 
vert  these  data  to  a  f^rm  suitable  for  actual  code  input. 

This  report  presents  a  detailed  illustration  of  the  analyses  involved 
in  deriving  the  soil  constitutive  properties  required  for  a  specific  code 
formulation  using  laboratory  test  data  from  just  one  stratum  of  a  single 
site.  Data  are  available  from  static  and  dynamic  uniaxial  strain  and  tri- 
axial  shear  tests;  code  property  requirements  are  for  mathematical  expres¬ 
sions  relating  mean  pressure  to  volumetric  strain,  Poisson's  ratio  to  mean 
pressure,  and  a  plastic  yield  criterion  to  mean  pressure. 

The  illustration  indicates  that  progress  is  being  made  in  developing 
mathematical  constitutive  models  that  are  realistic  in  terms  of  actual 
physical  behavior,  but  that  if  many  of  the  assumptions  presently  being 
made  in  soil  property  analyses  are  to  be  eliminated,  additional  soil  tests 
and  measurements  must  be  developed.  The  ill  '.stration  also  raises  questions 
as  to  the  validity  of  models  based  on  a  constant  Poisson's  ratio  or  a  con¬ 
stant  shear  modulus  and  suggests  that  the  behavior  of  the  various  models  be 
carefully  examined  under  different  states  and  paths  of  stress. 


ANALYSIS  OF  LABORATORY  TEST  DATA  TO  DERIVE 


SOIL  CONSTITUTIVE  PROPERTIES 


Introduction 


1.  The  primary  work  of  the  U.  S.  Army  Engineer  Waterways  Experiment 
Station  (WES)  Soil  Dynamics  Branch  for  the  Defense  Atomic  Support  Agency 
(DASA)  has  been  in  developing  soil  testing  facilities  and  in  conducting 
experimental  and  analytical  investigations  aimed  at  improving  both  the 
input  and  the  mathematical  formulation  of  constitutive  relations  used  in 
free-field  ground  shock  calculations .  This  is  an  area  of  considerable  cur¬ 
rent  interest  to  many,  as  numerous  site  investigations  that  require  the 
most  up-to-date  constitutive  property  definitions  and  computational  tech¬ 
niques  are  now  in  progress. 

2.  Last  year  (1968)  laboratory  test  data  from  several  sites  were 
presented  to  qualitatively  illustrate  the  effects  of  factors  such  as  load¬ 
ing  rate,  stress  history,  degree  of  saturation,  weathering,  and  state-of- 

stress  on  the  stress-strain  and  strength  properties  generally  used  in  vari- 

12 

ous  constitutive  property  formulations.  3  The  analyses  involved  in  de¬ 
riving  the  soil  constitutive  properties  required  for  a  specific  code  formu¬ 
lation  will  be  quantitatively  illustrated  herein  with  laboratory  test  data 
taken  from  just  one  stratum  of  a  single  site.  It  is  hoped  that  this  will 
permit  a  better  appreciation  of  the  uncertainties  and  assumptions  involved 
in  such  analyses  and  some  insight  into  the  research  work  ttiat  remains  to  be 
done . 

3.  Before  getting  into  details,  it  should  be  pointed  out  that  the 
particular  laboratory  tests  and  constitutive  model  that  will  be  used  in 
the  illustration  are  not  necessarily  the  most  recent  innovations.  There 
are,  as  would  be  expected  in  an  active  research  program,  other  tests  and 
other  constitutive  models  in  various  stages  of  development  and  evaluation. 
But  the  tests  and  model  that  will  be  used  in  the  illustration  presented 
herein  do  represent  the  current  state-of-the-art,  or  what  is  now  available 
for  immediate  application. 
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Basic  Constitutive  Relation  and  'Typical  Property  Requirements 


4.  Hie  basic  constitutive  relation  used  to  define  stress-strain  be¬ 
havior  utilizes  the  classical  linear  elastic  relation  between  the  stress 
tensor  and  the  strain  tensor  and  a  yield  criterion 

a.  .  =  Ke6.  .  +  2G  (e.  .  -■|e6..'\  and  .[jj.  =  f(p) 

U  v  10  3  10/  V  2max 

=  total  stress  tensor  =  p6 . .  +  a . 

ij  lj 

=  total  strain  tensor  =  ^  e6.  .  +  e' . 

3  10  iJ 

=  Kronecker  delta  function 

=  deviator  stress  tensor 

=  deviator  strain  tensor 

=  mean  normal  stress  =  ^ 

=  volumetric  strain  = 

=  bulk  modulus  -  ^ 

de  a:  . 

=  shear  modulus  =  g 

ij  i 

=  second  invariant  of  stress  deviation  =  —  o'.  .  o'.  . 

2  io  10 

The  constitutive  model  can  be  nonlinearized  by  using  incremental  stress- 
strain  relations  and  defining  one  or  both  of  the  soil  property  coefficients 
K  and  G  as  functions  of  one  of  the  stress  or  strain  invariants  such  as 
p  or  e  .  Inelasticity  can  be  incorporated  by  programming  two  sets  of 
property  coefficient  functions,  one  for  use  during  virgin  loading  and  one 
for  use  during  unloading  or  reloading. 

5.  Soil  property  input  to  a  computer  code  employing  this  basic  con¬ 
stitutive  model  can  be  specified  in  a  variety  of  forms.  One  typical  set  of 
constitutive  property  requirements  consists  of  (a)  hydrostat  expressions 
relating  mean  pressure  and  volumetric  strain  for  both  loading  and  unloading¬ 
reloading,  (b)  a  companion  set  of  expressions  for  Poisson's  ratio  v  as 


where 


ij 
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a  function  of  pressure,  and  (c)  a  yield  criterion  also  as  a  function  of 
pressure,  or 


PL  =  f^e)  and  Pu  =  ^(e.e^) 

vL  =  f3(p)  and  =  f4(p,Pinax)  where  v 


=  f5(p) 

6.  The  unloading-reloading  expressions  must  actually  define  families 
of  equations  that  are  dependent  on  previous  stress  or  strain  history.  To 
generate  these  families,  the  current  codes  accept  as  input  an  expression 
for  a  single  unl.oading-reloading  curve  and  then  apply  a  single-axis  trans¬ 
lation  to  the  various  points  of  unloading,  specified  as  either  e  or 

max 

p  .To  prevent  the  model  from  developing  energy-generating  hysteresis 
max 

loops  due  to  this  type  translation,  it  must  be  required  that,  at  any  given 
pressure,  the  unloading  bulk  modulus  be  equal  to  or  greater  tnan  the  virgin 
loading  modulus,  or 


de  —  de 


Equations  and  Assumptions  Used  for  Analyses  of  Laboratory  Tests 


7.  In  order  that  meaningful  soil  property  definitions  of  the  type 
just  described  can  be  obtained  from  the  results  of  laboratory  tests,  the 
tests  must  be  conducted  under  rigidly  controlled  states  of  stress  for  which 
all  components  of  the  stress  and  strain  tensors  can  be  defined,  by  measure¬ 
ment,  by  imposed  boundary  conditions,  or  by  assumption.  Before  the  various 
applied  loadings,  boundary  conditions,  and  measurements  applicable  to  the 
specific  tests  conducted  in  the  WES  laboratory  are  outlined,  the  general 
equations  and  assumptions  used  in  the  analyses  of  the  data  from  all  of  the 
laboratory  tests  should  be  examined. 
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Stress  and  strain  components 


8.  Total  stresses  within  the  soil  mass  beneath  the  site  being  in¬ 
vestigated  are  composed  of  initial  geostatic  or  overburden  stresses  and 
live  stresses  induced  by  the  blast  loading,  or 


^ TOTAL  =  0 OVERBURDEN  +  °LIVE  where  °TOTAL  -  0 

It  is  generally  assumed  that  the  soil  mass  cannot  support  tension  stresses, 
as  indicated  by  the  condition  that  the  total  stresses  must  al\vays  be  equ.i 
to  or  greater  than  zero.  Note  that  this  does  not  mean  that  live  tension 
stresses  cannot  be  supported;  they  can  be  supported  up  to  the  magnitudes 
of  the  overburden  stresses,  which  increase  continuously  with  depth.  Rea¬ 
sonable  estimates  can  be  made  for  the  magnitudes  of  these  overburden 
stresses,  which  permits  the  setting  of  tension  limits  for  the  live  stress 
components.  Since  the  overburden  stresses  are  relieved  when  the  soil  sam¬ 
ple  is  removed  from  the  ground,  the  test  specimens  must  be  statically 
recompressed  under  the  estimated  overburden  pressure  prior  to  application 
of  the  simulated  live  loading. 

9.  Total  strains  also  consist  of  overburden  and  live  components. 
But,  unlike  the  case  with  stresses,  overburden  strains  eQ  cannot  be 
readily  estimated  and,  for  all  practical  purposes,  are  indeterminate,  or 


eTOTAL  eOVERBURDEN  +  ®LIVE 


where  is  indeterminate 


Therefore,  all  WES  laboratory  stress- strain  and  strength  results  and  the 
constitutive  relations  derived  from  them  are  given  in  terms  of  live  stress 
and  live  strain.  This  means  that  the  zero  position  on  the  axes  of  the  var¬ 
ious  plots  is  assumed  to  represent  the  position  of  the  specimen  in  situ, 
and  the  plots  themselves  depict  response  of  the  specimen  to  applied  live 
loadings  or  departures  from  the  original  in  situ  stress  condition. 

General  equations  for  axial  symmetry 

10.  The  laboratory  tests  currently  used  for  WES  soil  property  in¬ 
vestigations  are  all  axially  symmetric  so  the  specimens  can  be  analyzed  in 


U 


terms  of  cylindrical  coordinates  using  the  two  general  equations  of 
equilibrium 


and  the  four  strain-displacement  relations 


__  du  _  u  _  Sw  .  _  1  /  du  ,  5w  \ 

er  ~  or  5  S6  ~  r  5  z  ~  bz  ’  rz  ~  2  V  Sz  Br  / 


provided  that  a  few  assumptions  are  made. 

Assumptions 

11.  First,  it  is  assumed  that  the  inertia  stresses  of  the  specimen 

52u  52w 

are  negligible,  which  allows  the  acceleration  terms  p  —r  and  p  — -  to 

be  dropped.  Insofar  as  possible,  this  condition  is  ensured  by  controlling 
the  loading  rise  times  in  the  experiments.  Next,  the  gravity  stresses  due 
to  the  weight  of  the  relatively  small  test  specimens  are  considered  negli¬ 
gible,  which  allows  the  y  term  to  be  dropped;  but,  as  was  pointed  out 
earlier,  the  total  gravity  or  overburden  stresses  acting  on  the  specimen 
are  not  negligible  and  must  be  applied  as  part  of  the  boundary  loading. 

12.  It  is  also  assumed  that  all  shear  stresses  and  all  shear  strains 
are  negligible,  which  permits  all  a  and  e  terms  to  be  dropped. 
Ensuring  this  condition  is  attempted  by  sealing  the  specimens  and  applying 
the  boundary  loading  through  fluids  wherever  possible.  Having  thus  elimi- 

ocr 

£ 

nated  all  but  the  term  in  the  second  equilibrium  equation,  it  must 

be  zero,  which  implies  that  a  uniform  state  of  vertical  stress  exists 
within  the  specimen. 

13.  Having  established  one  condition  of  uniform  stress,  it  is  simply 
assumed  that  all  other  states  of  stress  and  all  states  of  strain  are 
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14.  The  laboratory  data  for  constitutive  property  analyses  are 
obtained  from  two  basic  tests ,  i .  e . ,  the  uniaxial  strain  and  the  triaxial 
shear  tests.  A  brief  outline  of  the  boundary  and  loading  conditions, 
measured  responses,  and  the  types  of  data  plots  generally  obtained  from 
these  tests  follows.  Details  regarding  equipment  and  test  techniques  are 
given  in  references  4,  5,  6,  and  7- 

Uniaxial  strain  test 

15.  In  the  uniaxial  strain  test  (see  fig.  l)  a  condition  of  zero 
radial  strain  is  imposed  as  a  boundary  restraint  on  the  specimen  while  a 
controlled  vertical  stress  is  applied  to  it.  The  response  of  the  specimen 
is  measured  in  terms  of  a  vertical  surface  displacement  AL  ,  which  is  con¬ 
verted  to  Lagrangian  vertical  strain  e  by  dividing  the  displacement 

z 

time  history  by  the  original  height  of  the  specimen.  Since  the  radial 
strain  is  zero,  the  vertical  strain  determined  from  the  uniaxial  strain 
test  also  defines  volumetric  strain  e  .  Instantaneous  or  Eulerian  volu¬ 
metric  strain  can  be  obtained  from 

Ay  ez 
e  V  =  1  -  e 

z 

Plots  of  a  as  a  function  of  e  are  the  primary  end  product;  the  slopes 
z  z 

of  these  curves  define  the  constrained  modulus  M  .  The  constrained 
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Fig.  1.  Data  available  from  uniaxial  strain  test 

modulus  is  related  to  the  bulk  modulus  K  and  Poisson's  ratio  v  needed 
for  the  constitutive  model;  but  does  not  define  either  explicitly. 

1.6.  When  the  specimen  is  unloaded,  all  stress  (including  both  the 
live  stress  and  the  overburden  prestress)  is  removed  dynamically.  This 
permits  measurement  of  the  response  of  the  specimen  to  live  tension 
stresses  up  to  the  limit  set  by  the  overburden  as  shown  in  fig.  1.  By  in¬ 
cluding  this  negative  portion  of  the  unloading  curve  in  the  constitutive 
property  formulation,  the  necessity  for  including  a  gravity  term  in  the 
equations  of  motion  for  the  code  is  considered  to  have  been  eliminated, 
which,  if  correct,  should  simplify  the  computational  scheme. 

17.  Although  its  measurement  is  difficult,  the  radial  stress  re¬ 
quired  to  maintain  a  condition  of  zero  radial  strain  can  often  be  measured 
during  static  uniaxial  strain  tests.  Such  a  measurement  then  permits  plot¬ 
ting  o  as  a  function  of  o  .  These  plots  are  most  useful  in  that  from 

A  X 

them,  Poisson's  ratio,  mean  pressure,  and  the  square  root  of  the  second  in¬ 
variant  of  stress  deviation  can  be  calculated  directly,  since  for  the 
condition  of  uniaxial  strain 
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Aa 


v  = 


A  a  +  Aa 


T>  = 


+  2  a 
r 


a  -  a 
z  r 


Triaxial  shear  tests 

18.  In  the  triaxial  shear  test  (see  fig.  2),  a  constant  all-around 

stress  a  is  first  imnos^l  on  the  cnecimen.  Then  while  this  stress  is 
c 

maintained  as  a  boundary  condition  in  the  radial  direction  (i.e.,  o 

\  AL  r 

-  constant;,  either  controlled  axial  deformation  rates  ^  or  con¬ 
trolled  axial  stresses  are  applied  until  the  specimen  fails  in  shear. 
The  response  of  the  specimen  is  measured,  either  as  a  vertical  strain  for 
the  controlled  stress  tests  or  as  an  axial  deviator  stress  (or  principal 
stress  difference)  for  the  controlled  deformation  tests.  Tests  can  be 
conducted  both  dynamically  and  statically;  the  maximum  radial  stress  cur¬ 
rently  used  with  dynamic  tests  is  ’00  psi,v  but  static  tests  can  be  con¬ 
ducted  with  radial  stresses  up  to  10,000  psi. 


Fig.  2.  Data  available  from  triaxial  shear  tests 


*  A  table  of  factors  for  converting  British  units  of  measurement  to  metric 
units  is  presented  on  page  vii. 
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19.  In  either  case;  the  primary  product  is  a  plot  of  principal 

stress  difference  versus  vertical  strain  as  shown  in  fig.  2.  Only  the 

stress  information  at  "failure"  or  (a  -  a  )  as  a  function  of  a  is 

z  *  max  1 

currently  used;  these  data  are  used  to  plot  an  envelope  of  maximum  prin¬ 
cipal.  stress  difference  versus  pressure.  Such  an  envelope  is  directly 


proportional  to  the  yield  function 


Derivation  of  Constitutive  Properties 


20.  It  is  significant  to  note  that  although  the  uniaxial  strain  and 
triaxial  shear  tests  yield  some  obviously  useful  stress-strain  and  strength 
property  information,  neither  of  them  directly  yields  the  specific  proper¬ 
ties  outlined  in  paragraph  5  as  required  input  for  the  constitutive  model. 
It  is,  therefore,  appropriate  now  to  go  through  the  steps  involved  in  de¬ 
riving,  from  laboratory  test  dat  .  of  the  type  just  described,  a  complete 
set  of  constitutive  properties  in  a  form  suitable  for  code  input. 

Selection  of  site  profile 

21.  One  of  the  first  steps  is  the  selection  of  a  typical  site  pro¬ 
file.  No  two  boring  logs  are  ever  exactly  alike;  but  since  the  ground 
shock  codes  are  presently  only  ‘-'■''-dimensional,  a  single  boring  profile 
has  to  be  constructed  that  is  assumed  to  be  most  representative  of  the  area 
to  be  included  in  the  calculation.  Of  course,  the  more  borings  available 
for  such  an  analysis,  the  better. 

22.  The  soil  profile  shown  in  fig.  3  came  from  a  site  near  Valley 
City,  North  Dakota.  As  can  be  seen  in  fig.  3?  the  profile  is  made  up  of 
graphic  symbols  and  word  descriptions  that  are  used  as  guides  to  divide 
the  site  into  layers  or  zones;  a  set  of  meaningful  constitutive  properties 
must  then  be  defined  for  each  zone  selected.  Selection  of  the  layers  obvi¬ 
ously  has  to  be  coordinated  with  the  ground  shock  calculator,  since  the 
number  of  layers  and  the  minimum  layer  thickness  are  functions  of  the  code, 
the  computer  available,  and  the  time  step  to  be  used. 

23.  For  example  purposes,  a  10-ft-thick  zone  from  a  depth  of  4 5  to 
55  ft  below  the  ground  surface  was  selected.  The  blowup  of  this  zone  in 
fig.  3  shows  the  location  of  various  laboratory  test  specimens  within  the 
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zone.  For  the  example 
analysis,  there  are  re¬ 
sults  from  three  dynamic 
uniaxial  strain  tests, 
•jhree  static  uniaxial 
strain  tests  vrilth  radial 
stress  measurements,  six 
dynamic  triaxial  shear 
tests,  and  nine  static 
triaxial  shear  tests. 

Selection  of 
representative  data 
for  each  type  test 

24.  The  next  step 
is  to  select  representa¬ 
tive  data  from  this 
layer  for  each  type  test. 
This  is  by  far  the  most 
important  step  and 
should  be  much  more  than 
a  simple  quantitative 
averaging  of  the  test 
results.  Every  possible  check  should  be  made  to  ensure  that  the  selections 
represent  logical  conclusions  that  are  consistent  with  other  available  data 
and  at  least  with  the  more  basic  principles  of  soil  mechanics . 

25.  Results  of  the  three  dynamic  uniaxial  strain  tests  are  shown  in 
fig.  4;  the  curve  shown  by  the  dashed  line  is  considered  to  be  the  most 
representative  dynamic  uniaxial  strain  response  for  the  entire  10-ft-thick 
layer.  Also  shown  are  the  average  composition  properties  of  water  content 
w  ,  dry  unit  weight  V,  ,  specific  gravity  of  soil  solids  G  ,  void  ratio 

G.  S 

e  ,  percent  saturation  S  ,  and  the  percentages  by  volume  of  air  V  , 

0. 

water  V  ,  and  solids  V  . 
w  s 

26.  The  individual  test  specimens  were,  on  the  average,  99*3  percent 
saturated  and  contained  0.3  percent  air  by  vo?ume.  These  properties  refer 
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Fig.  3.  Soil  profile  and  test 
specimens  used  in  analyses 
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to  the  composition  of 
the  specimens  after  the 
static  overburden  stress 
has  been  reapplied; 
urior  to  this  loading 
the  specimens  in  the 
laboratory  were,  on  the 
average,  only  97.3  per¬ 
cent  saturated.  This 
2.0  percent  increase  in 
saturation  can  be  ac¬ 
counted  for  by  compres¬ 
sion  and  solution  of 
free  air  caused  by 
larger  pore  water 
stresses  in  the  in  situ 
condition  than  those  ex¬ 
isting  in  the  specimen 
upon  its  removal  from  the 
ground.  The  relation  be¬ 
tween  in  situ  degree  of 
saturation  and  its  value 
at  the  surface  is  given 

by  the  following 

8 

equation: 


AVERAGE  COMPOSITION 
PROPERTIES 
W  =23.1% 

—  0'd=  I  01.4  PCF  — 
GS  =  2.6I 

e  =  0.806 

S  =99.3%  ~J~ 


0.4  0.6 

% 


Fig.  4.  Dynamic  uniaxial  strain  test  results 


where 


1+  — 


1  +  i  (1  -  I,)So 
a 


S  -  in  situ  degree  of  saturation  (or  saturation  as  determined  from 
laboratory  testa  after  application  of  overburden  simulating 
preload  =  O.993) 

^o  “  of  saturation  at  ground  surface  (or  saturation  as  deter¬ 

mined  from  laboratory  tests  prior  to  application  of  overburden 
simulating  preload  =  O.973) 

u  =  in  situ  pore  water  pressure  (gage  pressure) 
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u  =  atmospheric  pressure  =  14.7  psi 

£L 

h  =  gas  solubility  constant  (0.02  for  air) 

Substitution  in  the  above  equation  shows  that  the  increased  saturation  is 
associated  with  a  pore  water  pressure  u  increase  of  approximately  13  psi; 
a  pore  pressure  increase  of  this  magnitude  would  require  that  the  specimens 
be  located  at  least  30  ft  below  the  free  water  table.  This  is  in  good 
agreement  with  the  field  data  and  suggests  that  the  specimen  compositions 
in  the  laboratory  prior  to  live  load  application  are  reasonably  represent¬ 
ative  of  those  of  the  in  situ  specimens. 

27.  Assuming  that  the  initial  live  strain  also  represents  air  being 

forced  into  solution,  100  percent  saturation  would  occur  at  a  live  vertical 

strain  of  0.3  percent  or  at  a  live  vertical  stress  of  200  psi,  as  shown  in 

fig.  4.  This  result  is  also  quite  reasonable  since  the  curves  show  a 

definite  stiffening  effect  at  that  point,  which  is  characteristic  of 
2  9 

saturation.  57 

28.  For  this  example,  the  representative  dynamic  uniaxial  strain 
relation  has  been  defined  to  a  maximum  vertical  stress  of  2000  psi.  This 
cannot  be  a  purely  arbitrary  decision  on  the  part  of  the  soils  engineer, 
but  must  be  worked  out  in  advance  with  the  calculator  based  on  the  maximum 
pressure  to  be  input  at  the  starting  ground  range  and  preliminary  estimates 
of  stress  attenuation  with  depth. 

29.  Results  of  the  three  static  uniaxial  strain  tests,  with  radial 

stress  as  well  as  vertical  strain  measurements,  are  given  in  fig.  5.  Data 

for  these  tests  were  available  up  to  a  vertical  stress  of  only  about 

350  psi.  The  most  representative  relations  between  a  and  e  and  be- 

z  z 

tween  a  and  cr  are  defined  by  the  dashed  lines.  According  to  the 
z  r 

specimen  composition  data  shown  in  fig.  4,  100  percent  saturation  should 
occur  at  a  strain  of  0.3  percent  or  at  a  vertical  stress  of  about  35  psi 
and  a  radial  stress  of  about  20  psi;  this  is  equivalent  to  a  mean  pressure 
of  only  25  psi.  For  the  dynamic  tests,  the  vertical  stress  required  for 
full  saturation  was  200  psi  or  an  estimated  mean  pressure  of  about  150  psi. 

30.  Fig.  6  shows  static  triaxial  shear  test  results  for  four 
controlled  rate  of  deformation  tests  with  values  up  to  70  psi  and  five 
controlled  stress  tests  with  values  up  to  1470  psi.  The  maximum 
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Fig.  5.  Static  uniaxial  strain  test  results 


Fig.  6.  Static  triaxial  shear  test  results 
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Fig.  7»  Dynamic  triaxial  shear  test  results 


principal  stress  differ¬ 
ences  from  these  tests  have 
been  plotted  as  a  function 
of  mean  pressure,  and  the 
most  representative  static 
failure  envelope  was  drawn 
for  the  layer.  Note  that 
the  envelope  is  flat,  or  of 
the  von  Mis  is  type,  through 
its  entire  pressure  range. 
This  indicates  that  satura¬ 
tion  occurred  at  a  very  low 
pressure,  which  is  in  agree¬ 
ment  with  the  static  uni¬ 
axial  strain  test  data 
sho>m  in  fig.  5. 

31.  Fig.  7  shows  the 
results  from  six  dynamic 
triaxial  shear  tests  con¬ 
ducted  with  radial  stresses 
up  to  only  170  psi.  The 
results  with  these  rela¬ 
tively  low  a  values  in- 
r 

dicate  a  gradual  increase 
in  maximum  principal  stress 


difference  with  mean  pressure,  but  the  most  representative  envelope  has 


been  flattened  at  a  pressure  of  about  150  psi  since  the  dynamic  uniaxial 
strain  test  results  indicated  100  percent  specimen  saturation  at  that  pres¬ 
sure.  The  maximum  dynamic  strength  is  62  psi,  which  is  I.55  times  the 


maximum  static  value  of  ho  psi  shown  in  fig.  6. 

Construction  of  property  relations 

?2.  After  the  representative  test  data  have  been  selected,  the 
final  step  is  construction  of  the  three  constitutive  property  relations  as 
outlined  for  computer  Input:  pressure  versus  volumetric  strain  functions, 


lh 


Poisson's  ratio  versus  pressure  functions,  and  a  yield  function. 

33.  The  first  to  be  determined  is  the  yield  function.  A  limiting 

principal  stress  difference  versus  pressure  envelope  from  the  static  tri-  - 

axial  test  results  was  defined  (see  fig.  6),  but  as  noted  in  paragraph  17, 

principal  stress  difference  versus  pressure  for  a  complete  load-unload 

cycle  can  be  calculated  directly  from  the  results  of  static  uniaxial  strain 

tests  with  radial  stress  measurements.  The  a  versus  a  values  from 

z  r 

the  plot  given  in  fig.  5  as  the  most  representative  static  uniaxial  strain 
response  for  the  layer  have  been  converted  to  values  of  (a^  "  ar)  versus 
p  in  table  1  so  that  a  comparison  can  be  made  with  the  static  triaxial 
test  results. 

34.  The  comparison  between  the  static  uniaxial  strain  test  relation 
for  (a  -  a  )  versus  p  and  the  static  triaxial  shear  test  results  is 
shown  in  fig.  8.  The  heavy  solid  line  at  a  principal  stress  difference 

of  40  psi  is  the  representative  triaxial  shear  envelope  determined  from  the 
static  test  dataj  its  mirror  image  is  plotted  at  -40  psi  as  a  limit  for 
unloading.  The  other  solid  line  is  the  path  of  principal  stress  difference 
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Fig.  8.  Static  uniaxial  strain  and  triaxial  shear’ 
results  adjusted  for  compatibility 
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versus  pressure  plotted  from  the  static  uniaxial  strain  results  given  in 
table  1.  Note  that  a  constant  principal  stress  difference  of  25  psi,  or 
yield,  is  reached  at  a  pressure  of  about  100  psi.  Yielding  continues  until 
the  maximum  test  pressure  is  reached;  then  the  path  abruptly  leaves  the 
upper  yield  surface  on  unloading,  but  appears  to  yield  again  along  a  lower 
yield  surface  after  the  pressure  decreases  to  about  200  psi. 

35.  Since  the  triaxial  test  data  indicated  a  yield  envelope  with  a 
maximum  principal  stress  uifi'erence  of  _^40  psi  and  the  uniaxial  strain  data 
gave  +25  psi,  some  adjustment  had  to  be  made  to  make  them  compatible.  More 
weight  was  given  to  the  trlaxial  test  data  in  this  case,  and  a  maximum 
yield  value  of  35  psi  war  selected,  .lie  heavy  dashed  line  in  fig.  8  repre¬ 
sents  the  downward  adjust  jd  triaxial  shear  envelope  and  the  lighter  dashed 
line,  the  upward  adjusted  path  of  uniaxial  strain.  Both  the  loading  yield 
pressure  value  of  100  psi,  at  point  3.,  and  the  unloading  value  of  200  psi, 
at  point  3,  were  retained.  The  pressure  at  point  2  is  determined  by  the 
maximum  vertical  stress  of  350  psi,  and  that  at  point  4  represents  the 
limit  for  unloading  set  by  the  weight  of  the  overburden  as  a  vertical 
stress  of  -50  psi. 

36.  The  static  uniaxial  strain  and  triaxial  shear  results  have  now 
been  adjusted  to  be  compai  iblo  with  each  other.  The  next  step  is  to  adjust 
the  dynamic  shear  envelope  to  be  compatible  with  the  adjusted  static  shear 
envelope,  'he  dashed  li.ie  L;i  tig.  9  is  the  static  shear  envelope  that  was 
adjusted  downward  from  a  nax'mum  prit cipal  stress  difference  of  40  psi  to 
35  psi.  The  dynamic  envelope,  shown  by  the  solid  line,  has  also  been 
adjusted  downward  from  a  test-determined  maximum  principal  stress  differ¬ 
ence  value  of  62  osi  tc  a  now  value  of  55  psi  in  order  to  preserve  the 
1.55  ratio  between  dyna'uic  yield  strength  and  static  yield  strength  indi¬ 
cated  by  the  laboratory  tests. 

37-  The  dynamic  .shear  envelorx  for  (o  -  o  )  as  a  function  of 

max 

p  has  now  seen  established  for  the  constitutive  model  and  can  be  readily 
fit  idth  a  polynomial  to  give  a  computer-acceptable  equation  for  the 
dynamic  yield  relation  as  shown  in  fig.  9- 

38.  Next,  a  Poisson's  ratio  versus  pressure  relation  is  required. 
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ADJUSTED  DYNAMIC  SHEAR  ENVELOPE 


MEAN  NORMAL  PRESSURE  p,  PSI 
EQUATION  FOR  DYNAMIC  YIELD  RELATION 


-/Ji  = - =  25.8  r  0.0953  p  - 0.000370p*  FOR  p<  150  PSI 

MAX  1 3 

=  31.8  FOR  pi  150  PSI 

Fig.  9.  Adjusted  dynamic  shear  envelope  and  yield  equation 


In  order  to  get  that,  a  path  of  principal  stress  difference  versus  pres¬ 
sure  for  dynamic  uniaxial  strain  must  be  established  to  go  with  the  ad¬ 
justed  dynamic  shear  or  yield  envelope.  Since  the  necessary  dynamic  data 
for  this  stress  path  are  not  available,  one  has  to  be  constricted  based  on 
static  information,  'the  heavy  solid  lines  in  fig.  10  represent  the 


mmmm 


dynamic  shear  envelope  and  the  lighter  solid  lines,  the  static  shear  en¬ 
velope.  The  light  dashed  line  shows  the  (a  -  a  )  versus  p  path  for 

Z  X* 

static  uniaxial  strain  during  loading  to  yield  and  unloading  between  yield 
surfaces.  The  slopes  along  this  path  are  related  to  Poisson's  ratio  by 
the  expression 

d(az  "  CTr ^  _  3(1  -  2 v) 
dp  -  1  +  v 

A  slope  cf  3  represents  a  Poisson's  ratio  of  zero,  and  a  zero  or  horizontal 
slope  represents  a  Poisson's  ratio  of  0.5. 

39*  Experience  has  been  that,  although  the  static  and  dynamic  soil 
test  responses  may  be  quantitatively  different,  they  are  quite  similar  in 
form  as  indicated  by  the  results  given  in  figs.  4,  5>  6,  and  7.  Thus,  the 
first  trial  dynamic  stress  path  of  (a  -  a^)  versus  p  (shown  by  the  heavy 
dashed  line  in  fig.  10  between  the  origin  and  point  1  and  between  points 
2  and  3)  was  given  a  shape  similar  to  the  static  path.  With  this  path, 
Poisson's  ratio  continuously  increases  over  a  pressure  range  of  150  psi 
(based  on  the  estimated  dynamic  uniaxial  strain  saturation  pressure)  during 
both  loading  and  unloading.  A  much  more  simplified  assumption  is  to  let 
Poisson's  ratio  be  constant.  A  second  trial  path  was  constructed  for  a 
constant  Poisson's  ratio  of  0.25  as  shown  by  the  dotted  lines  in  fig.  10. 

40.  Wow,  some  assumption  must  be  made  about  unloading  from  maximum 
pressures  other  than  that  associated  with  a  vertical  stress  cf  350  psi. 

The  most  logical  guess  would  be  a  simple  p-axis  translation  as  indicated  by 
the  heavy  dashed  lines  in  fig.  10  constructed  parallel  to  the  one  between 
points  2  and  35  such  a  translation  is  consistent  with  a  constant  Poisson's 
ratio  or  single-slope  assumption. 

^1.  Once  two  trial  dynamic  stress  paths  of  principal  stress  differ¬ 
ence  versus  pressure  for  a  state  of  uniaxial  strain  have  been  established, 

they  can  readily  be  converted  to  more  useful  plots  of  a  versus  a  . 

z  r 

The  long-dashed  line  in  fig.  11  shows  the  adjusted  static  uniaxial  strain 
relation  between  and  ar  .  The  solid  line  shows  the  first  trial  dy¬ 
namic  relation:  point  1  indicates  the  point  of  first  yield  during  loading; 
point  2  corresponds  to  the  maximum  vertical  stress  of  350  psi;  point  3 


represents  first  yield 
during  unloading;  and 
point  4  the  limiting  ten¬ 
sion  value  of  ct  =  -50 
z  y 

psi.  The  data  for  this 
trial  dynamic  stress  path 
are  given  in  table  2. 

The  short-dashed  lines  in 
fig.  11  show  the  second 
trial  or  idealized  dy¬ 
namic  uniaxial  strain 
relation  corresponding 
to  a  constant  Poisson's 
ratio  of  0.25. 

42.  Nov/  that  there 
are  dynamic  uniaxial 

strain  values  for  o 

z 

versus  a r  as  well  as 
for  az  versus  e  ,  a 
curve  of  mean  pressure 
p  versus  volumetric 


Fig.  11.  Vertical  stress  versus  radial 
stress  for  uniaxial  strain  based  on 
trial  stress  paths 


strain  e  can  be  plotted.  But  since  the  plot  must  extend  to  a  vertical 
stress  of  2000  psi,  the  az  versus  or  relation  to  this  maximum  stress 
must  first  be  extended  according  to  the  p-axis  translation  assumption.  The 
tabulated  data  for  this  extended  az  versus  a,  relation  from  the  first 
trial  dynamic  uniaxial  strain  stress  path  and  the  corresponding  p  versus 
e  function  are  given  in  table  3.  The  light  solid  line  in  fig.  12  shows 
the  az  versus  a.  relation,  first  unloaded  from  a  a,  of  350  psi  and 
then  extended  to  a  a2  of  2000  psi  and  unloaded.  The  corresponding  p 
versus  e  relation  is  shown  by  the  heavy  solid  line. 


43.  This  trial  combination  of  az  and  a,  produces  an  interesting 
result  in  the  p  versus  e  relation  at  high  pressure;  the  unloading  vol¬ 
umetric  strains  are  less  than  the  loading  strains.  This  portion  of  the 
data  has  been  plotted  on  an  expanded  scale  in  the  upper  left-hand 
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Fig.  12.  Extended  oz  versus  ar 
p  versus  e  based  on  trial  stress 
paths  for  uniaxial  strain 


corner  of  fig.  12.  Again,  the 
light  solid  line  is  the 
versus  a  relation,  and  the 
heavy  solid  line  is  the  p 
versus  e  relation  shown  loop¬ 
ing  backwards  as  it  unloads. 

The  dashed  lines  were  plotted 
from  data  taken  from  the  trial 
in  which  Poisson's  ratio  was 
equal  to  0.25  and  yield  essen¬ 
tially  the  same  p  versus  e 
curve  as  the  first  trial  data. 

44.  Since  a  model  that 
generates  strain  energy  instead 
of  absorbing  it  was  not  con¬ 
sidered  to  be  a  good  solution, 
it  was  decided  that  the  exten¬ 
sion  of  the  a  versus  a 
z  r 

relation  to  higher  pressures  " 
might  have  been  done  a  little 
hastily.  Therefore,  determina¬ 
tion  of  the  o  versus  o 
z  r 

relation  that  would  eliminate 
this  strain  reversal  in  the 


p  versus  e  function  at  high 

pressure  by  permitting  the  p  versus  e  curve  to  at  least  initially  un¬ 
load  along  its  original  loading  path  was  necessary.  Such  a  relation  can  be 
calculated  from  a  tabulation  of  the  a  versus  e  values  from  the  repre- 

sentative  dynamic  uniaxial  strain  relation  given  in  fig.  4  and  the  corre¬ 
sponding  p  values  from  fig.  12  required  to  eliminate  the  reverse  hyster¬ 
esis  .loop  in  the  p  versus  e  relation.  Results  given  in  table  4  and 
fig.  13  show  that  when  unloading  takes  place  after  an  extended  loading  along 
a  yield  surface,  Poisson's  ratio  may  not  drop  abruptly  to  a  relatively  low 
/alue  such  as  0.25,  but  may  gradually  decrease  starting  with  a  much  higher 


20 


value  such  as  0.47  or  0.1*8. 

1*5.  Naturally,  the 
first  step  was  a  review  of 
the  pertinent  literature  to 
determine  if  some  physical 
evidence  could  be  found  to 
support  the  above  hypothe¬ 
sis.  An  excellent  example 
was  found  in  the  results  of 
a  static  uniaxial  strain 
test  with  radial  stress 
measurements  conducted  on 
Suffield  silty  clay  at  the 
University  of  Illinois;10 
these  test  results  are 
plotted  in  the  insert  at 
the  upper  left-hand 
corner  of  fig.  13.  The 
data  show  loading  to  a 
vertical  stress  of 
300  psi,  followed  by  an 
unloading  cycle  that  starts 
with  a  relatively  low  Pois- 


Fig.  13.  Adjusted  dynamic  vertical  stress 
versus  radial  stress  for  uniaxial  strain 


son  s  ratio  as  did  the  WES  static  test  data  and  a  subsequent  second  loading 
cycle  along  an  apparent  yield  surface  to  a  much  higher  vertical  stress, 
from  which  it  unloads  with  a  relatively  high  Poisson's  ratio,  i.e.,  between 
0.47  and  0.48. 


46.  With  this  bit  of  physical  evidence  as  backup,  the  adjusted 
dynamic  uniaxial  strain  relation  between  0^  and  <7  given  in  fig.  13  was 
used  in  the  constitutive  property  calculations.  This  plot  defines  Poisson' 
ratio  only  for  loading  from  the  origin  and  unloading  from  a  values  of 
350  and  2000  psi;  the  next  step  is  to  define  a  general  relation  for 
Poisson  s  ratio  that  accounts  for  unloading  from  any  maximum  pressure  apt 
to  be  encountered  during  the  ground  shock  calculation.  To  assist  in  this 


VIRGIN  LOADING 

•0=0.25  FOR  ALL  VALUES  OF  p 
UNLOADING  -  RELOADING 

FOR  0<pMAX<  1230  PSI 
■0=0.250 

FOR  Pmax>I230  PSI 
■0=0.476  WHEN  p>l340-0.0685  pMAX 
■0=0.250  WHEN  p£l340-0.0685  pMAX 
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Fig.  14.  Adjusted  and  idealized  stress  paths  for  uniaxial  strain 

step,  the  a  versus  a  relation  in  fig.  13  was  converted  to  the 

(az  -  ar)  versus  p  relation  she  .1  in  fig.  14.  The  heavy  lines  in  the 
upper  figure  show  the  loading  path  from  the  origin  and  the  unloading  path 
f0r  a  Pmax  value  corresponding  to  a  a  of  2000  psi  followed  by  unload¬ 
ing  along  the  recently  adjusted  path. 

47 .  Definition  of  the  actual  functional  form  for  these  stress  paths 
between  yield  surfaces  and  their  obvious  transitions  from  one  shape  to 
another  as  pressure  increases  is  an  item  that  requires  considerably  more 
testing  and  study.  For  the  present,  Poisson's  ratio  has  simply  been  de¬ 
fined  as  having  two  values,  0.250  and  0.476,  as  indicated  by  the  dashed 
lines  in  fig.  14.  The  value  of  0.250  is  used  during  any  virgin  loading  and 
WhGn  pmax  is  iess  than  123°  Psi5  for  Pmax  values  above  1230  psi,  the 
Poisson's  ratio  value  of  0.470  is  used  initially  and  is  changed  to  the 
0.2')0  value  when  the  unloading  path  defined  by 

p  =  p  -  564  +  10.250(o  -  o  ) 

max  z  r 
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intersects  the  boundary  defined  by 


p  =  1185  +  0.833 (ffz  -  or) 

Equations  can  be  -written  now  for  dynamic  Poisson's  ratio  as  a  function  of 
mean  pressure,  as  shown  in  the  upper  right-hand  corner  of  fig.  14. 

48.  Once  the  dynamic  yield  envelope  and  Poisson's  ratio  functions 
have  been  established,  they  can  be  combined  with  the  representative  dynamic 
ct  versus  e  relation  to  plot  the  final  dynamic  p  versus  e  relation 

Z  2 

shown  by  the  heavy  solid  line  in  fig.  15.  However,  the  unloading  slopes 
for  this  relation  at  pressures  below  about  100  psi  are  less  than  those 
during  virgin  loading,  which  violates  one  of  the  conditions  for  the  model 
given  in  paragraph  6.  To  avoid  this  violation,  the  unloading  "tail"  has 
been  stiffened  slightly  as  shorn  by  the  dashed  line.  The  final  step  is 
to  fit  both  loading  and  unloading  curves  with  polynomial  equations  as  shown 
at  the  top  of  fig.  15 . 


Ixg.  15.  Adjusted  dy¬ 
namic  pressure  versus 
volumetric  strain 
relation 
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49-  Finally,  the  completed  constitutive  model  is  in  equation  form: 


for  virgin  loading 

4  7  ?  10  ^ 

(1)  pT  =  6.1169  x  10  e  -  3-0821  X  10 ‘e  +  1.2467  X  10  e° 

L 

-  1.1060  X  lO^e4  +  3.2134  X  1013e5 

(2)  v_  =  0.250  for  all  p 

L 

(3)  */2  =  25.8  +  95-3  X  10“ 3p  -  37-0  X  10"5p2  for  p  <  150  psi 

V  max 

=31.8  psi  for  p  >  150  psi 

for  unloading-reloading 

(1)  py  =  7-0249  x  104(e  -  0.0015)  +  7.7503  x  106(e  -  0.0015)2 

+  4.2521  X  109(e  -  0.0015 )3 

-  I.6526  x  10i;L(e  -  0.0015 )4 

-  1.8721  X  1013(e  -  0.0015)5 

(2)  v  =  0.250  for  0  <  p  <  1230 

u  max 

=  0.476  for  p  >  1230  psi  and  when  p  >  1340  -  0.0885  p 

max  max 

=  0.250  for  p  >  1230  psi  and  when  p  <  1340  -  O.O885  p 

max  max 

(3)  /jI  =  -25.8  -  95-3  X  10" 3p  +  37-0  x  10-y  for  p  <  150  psi 

V  max 

=  -31-8  psi  for  p  >  150  psi 


Conclusion 


50.  At  the  beginning,  it  was  stated  that  a  better  appreciation  for 
the  uncertainties  in  a  constitutive  property  analysis  might  be  gained  from 
the  study  reported  herein.  The  many  assumptions  and  adjustments  show  that 
the  property  analysis  does  not  lead  to  a  "closed-form  solution."  BUt  an 
attempt  was  made  to  keep  things  logical  and  orderly  and  within  the  realm  of 
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physical  evidence  and  reasoning,  and  that  in  itself  is  certainly  a  step  in 
the  right  direction.  It  should  now  be  clear  that  there  is  considerably 
more  involved  in  formulating  constitutive  properties  than  just  conducting 
a  few  tests  and  plotting  the  results.  It  should  also  be  obvious  that  there 
must  be  a  high  degree  of  cooperation  and  interchange  between  the  soils 
engineers  and  the  calculators,  and  that  there  are  a  number  of  areas  of 
research  still  open  to  both. 

51.  For  instance,  as  far  as  laboratory  tests  are  concerned,  more 
strain  data  are  needed.  Only  one  strain  measurement  was  used  during  the 
entire  example  analysis :  that  was  the  vertical  strain  from  the  uniaxial 
strain  tests.  Measuring  soil  strain  is  very  tedious  and  difficult,  but 
progress  is  being  made  in  obtaining  radial  strain  measurements  during  tri- 
axial  compression  tests  which  promises  to  make  possible  the  measurement  of 
several  useiul  modulus  properties. 

52.  More  dynamic,  high-pressure  modulus  data  are  also  needed;  for 
while  measurement  of  constrained  modulus  has  been,  and  most  likely  will 
continue  to  be,  of  primary  importance,  other  coefficient  relations  such  as 
shear  modulus  and  Poisson’s  ratio  are  equally  important  for  general-purpose 
two-dimensional  models . 

53 .  This  leads  to  a  conclusion  for  the  calculators  that  the  constant 
v  and  constant  G  type  models  should  be  examined  carefully  in  order  that 
behavior  under  different  states  and  paths  of  stress  can  be  observed.  The 
unloading-reloading  logic  used  in  the  code  should  be  inspected  to  determine 
specif i cal] v  what  type  translations  are  made  to  accommodate  unloading  from 
any  maximum  signal  and  the  method  used  for  treating  reloading.  Another 
point  that  needs  investigating  is  whether,  in  fact,  gravity  stresses  can  be 
accounted  for  in  the  calculations  by  including  their  effects  in  the  consti¬ 
tutive  model  rather  than  by  including  an  added  term  in  the  equations  of 
motion . 

54.  It  is  hoped  that  these  suggestions  for  additional  tests  and  for 
perhaps  some  improvements  in  the  constitutive  models  have  not  left  the  im¬ 
pression  that  simply  having  more  test  data  of  more  different  types  and  hav¬ 
ing  more  generalized  constitutive  models  will  simplify  the  job  of  defining 
constitutive  properties  for  code  input.  Just  the  opposite  will  probably  be 
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the  case;  i.e.,  the  job  will  be  much  more  complicated  than  it  is  already. 
This  should  not  deter  such  efforts,  for  if  the  soils  engineers  and  the  cal 
culators  will  continue  to  work  together,  it  will  lead  to  more  realistic 
models  based  on  more  physical  facts,  and  that  cannot  help  but  lead  to 
better  ground  shock  calculations. 
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